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6.5 Geothermal Plume 

6-5g-plume-L.tex 

R..A.Wooding, (1963), J Fluid Mech. 15, 527-544. 

C. S. Yih, (1965), Dynamics of Nonhomogeneous Fluids, Macmillan. 

D. A. Nield and A. Bejan, (1992), Convection in Porous Media. Springer- Verlag. 
Consider a steady, two dimensional plume due to a source of intense heat in a porous 

medium. From Darcy's law: 

%»-- 1 < MJ > 

where k denotes the permeability, and 

V dp tar o\ 

- W = ---pg (6.5.2) 

These are the momentum equations for slow motion in porous medium. Mass conservation 
requires 

u x + w z = (6.5.3) 

Energy conservation requires 

dT OT fd 2 T d 2 T\ 

where 



U dx ~*~ W dz a \ dx 2 dz 2 ) (6.5.4) 



a = ^ (6.5.5) 

denotes the thermal difusivity. 
Equation of state: 

p = po (1 - 0(T - T )) (6.5.6) 
Consider th flow induced by a strong heat source. Let 

T-T = T, p = Po + p' 

where p is the hydrostatic pressure satisfying 

dpo n 
— « Pod = 0. 

Eqn. (6.5.2) can be written 

k W = ~% +9Po(3T '- (6 - 5 ' 7) 
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6.5.1 Boundary layer approximation 

Eliminating p' from Eqns. (6.5.7) and (6.5.1), we get 

^(w x -u z ) = gpo(3T' x . 

Let tp be the stream funciton such that 

u = if> z , w = -if) x 

then 

i> x:c + A z = -—T> (6.5.8) 
A* 

For an intense heat source, we expect the plume to be narrow and tall. Let us apply the 
boundary layer approximation and check its realm of validity later, 

d d 
u<^w, — > — . 

ox dz 

hence 

Yxx — ± x 

A* 

or 

^ s _£Po^ T ' ; (6 .5. 9) 
A 4 

which is the same as ignoring dp' /dz in Eqn. (6.5.7). 

This can be confirmed since u <C w dp' /dx ~ 0, p' inside the plume is the same as that 
outside the plume. But 

* 

oz 

outside the plume, hence dp' /dz inside as well. 
Applying the B.L. approximation to Eqn. (6.5.4) 

uT' x + wT' z = aT' xx (6.5.10) 

Using the continuity equation we get 

(uT% + (wT% = aT' xx . 

Integrating across the plume, 

d f°° 

— wT'dx = (6.5.11) 
since T' — outside the plume. It follows that 

/OO /"OO 
wT'dx = -p C if) x T'dx = Q = constant. (6.5.12) 
OO J — OO 
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6.5.2 Normalization 

Let us take 

x = Bx, z = Hz, u = ^-u, w = Ww, T -> AT6 (6.5.13) 

H 

where H, B, AT and W are to be determined to get maximun simplicity. We then get from 
the momentum equation, 



from the energy equation, 



_ T gpo(3AT 



aH 

u6 s + w9 z = yy^Oxx, 



and from the total flux condition, 

p CWBA j icOtl.r - () 

Let us choose 



9 "» flAT .1 (6.5.14) 



/j,W 
aH 



= 1 (6.5.15) 



WB 2 
and 

PoCWBAT = Q, (6.5.16) 
which gives three relations among four scales, B,H,W, AT. Then 

w = ijj x = -e, (6.5.17) 

from the energy equation, 

u0 £ + w6 E = 9 Sx , (6.5.18) 

and from the total flux condition, 

/oo 

w9dx = l (6.5.19) 
-oo 

In addition we require that 

w(±oo,2) = 0, 6(±oo,z)=0 (6.5.20) 

u(0,z) = dw ^ z ^ =0) x = . (6.5.21) 
ox 

From here on we omit overhead bars in all dimensionless equations for brevity. 
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6.5.3 Similarity solution 

Now let 

x = x a x * z = x b z* i) = \ c %i)* e = x d e*. 

FromEqn. (6.5.17) 



A c-a (W_\ = _ xd0 * 



For invar iance we require, 
From (6.5.19) 

therefore, 

From Eqn. (6.5.18) 

implying, 

Finally 



/ 



dx* 

c-a = d. (6.5.22) 

^-dx* x c - a+a+d = l. 
ox* 

a + d = 0. (6.5.23) 



yc+d—a—b y^d- 



2a 



c + a -b = 0. (6.5.24) 



a , a , 3 
c = -, a = —-, b = -a. 

2' 2' 2 



In view of these we introduce the following similarity variables, 

^ = ^2/3' i> = z 1/3 f(v), e = z- l ^h{ri). (6.5.25) 
Note that at the center line rj = 

w = -ip x oc 2 1 / 3 /'(0)(-)^ 2/3 ~ 2- 1 / 3 /'(0) ~ ^" 1/3 (6.5.26) 

oc z" 1/3 /i(0) (6.5.27) 

and 

6 oc z 2/3 (6.5.28) 

Thus the velocity and temperature along the centerline decay as z -1 / 3 and the plume width 
grows as z 2 / 3 . 

Substituting these into Eqns. (6.5.17) and (6.5.18), we get, after some algebra 



~ = h (6.5.29) 
drj 



and 



Tn Uh) = 3^. (6.5.30) 
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The boundary conditions are, 

/ = o (v = o) 

/"(0) = 0, (w(0, z) = w max ) 
/(±oo), /'(±oo) = 
h(±oo) = 0. 

Integrating Eqn. (6.5.30), we get 

fh = 3ti. 

Using Eqn. (6.5.29), we get 

//' = 3/". 

Integrating again, we get 

-6/ = / 2 - f 
where /„ = /max- Let / = -f F, then 

/ (l-F 2 ) = 6F',or 

which can be integrated to give 



2x art dF fo d7 l 



1 - F 2 6 



/o?7 = 1 1 + F 
6 2 n i-F 



Thus 



1 + F XV2 



or 



1 - F 

1 + F 
1 - F 



e hv/e 



g/orj/3 



Solving for F, we get 



e kl'i + i 6 
What is /„? Let us use Eqn. (6.5.29) 



e /o/3 _ i f „ 

-tanh^ (6.5.31) 



/.oo fir/ ./-oc 



since 

/' = -/oF' = -4sech 2 ^ 
6 6 

and 

h = -/'. 
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Therefore, 

f2\ 2 



Since 



we get /o! 



The solution is 



/oo 
sech 4 zdz = 4/3. 
■oo 



and 



'9\V3 

f ( - J ( 6 - 5 ' 32 ) 



1/3 /qX 1 / 3 
/Hi) tanh - | (6.5.33) 



& = -/' = -( II sech 2 f ^ 13 ^ (6.5.34) 



2) " \2) 6 



Computed results are given in Figures. 

Remark Checking the boundary layer approximation. 

d 2 ib 1 d 2 ib 
5a; 2 ' dz 2 

dx 2 dz 2 
hence for large z, B. L. approximation is good. 



V = i (6-5.35) 



6.5.4 Return to physcial coordinates 

Start from 

" " z 2 /3 

z 1 ^ = (6.5.37) 

Then 



x/£ /tf 2 / 3 \ / x 

?7 = 



(z/F) 2 / 3 V B 

By eliminating and AT from(6.5.35) and (6.5.37), we get 



) (M (^-38) 



CB 
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Figurs ft.fi. f; Theoretical solution for a geotfiermai pfume 

due to Yin. 

(Adapted from Yih, Dynamics of Nanhomogerteous Fiuids r 1965). 



From (6.5.36), we get 



It follows that 



Now 



H W 1 Qg/3 



B 2 a aV CB 



H 1 Qg/3 



B 3 / 2 a\ C 



fV3 VFS Vi?7 V WB I Kz 1 / 3 



(6.5.39) 



(6.5.40) 
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It can be shown that 

#1/3 / 1/3 ^ / g \ 1/3 

WB~ [q^_\B^) ~a^\Q^3J 



which depends on the fluid properties and the given heat source strength. 
Also 

*V30 = h ( v ) = [Hz) 1,z ATT' = (H^ 3 AT)z 1/3 T' (6.5.41) 

We can show that 

1/3 



which also depends on the fluid properties and the given heat source strength. 



